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Abstract In this paper we identify various inaccuracies in the paper by R. R. 
Saxena and S. R. Arora, A Linearization technique for solving the Quadratic 
Set Covering Problem, Optimization, 39 (1997) 33-42. In particular, we ob¬ 
serve that their algorithm does not guarantee optimality, contrary to what is 
claimed. Experimental analysis has been carried out to assess the value of this 
algorithm as a heuristic. The results disclose that for some classes of problems 
the Saxena-Arora algorithm is effective in achieving good quality solutions 
while for some other classes of problems, its performance is poor. We also 
discuss similar inaccuracies in another related paper. 

Keywords 0-1 programming, quadratic programming, algorithms, set 
covering problem, and heuristic. 


1 Introduction 

The set covering problem is well studied in the Operations Research literature 
miiiiiiiiiia. Most of the works on the problem reported in the literature have 
a linear objective function. Bazaraa and Goode introduced the quadratic 
set covering problem (QSP) and proposed a cutting plane algorithm to solve 
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it. Adams [T] and Liberti m proposed linearization techniques for binary 
quadratic programs. Since QSP is a binary quadratic programming problem, 
these linearization techniques can be used to formulate QSP as a 0-1 integer 
linear program. QSP is known to be NP-hard and polynomial time approxi¬ 
mation algorithms are also available [7] to solve special classes of this problem. 

Saxena and Arora m studied QSP and discussed various structural prop¬ 
erties of the problem along with a linearization algorithm which is claimed to 
produce an optimal solution. The notion of linearization used in HZ] is differ¬ 
ent from the concept of ” linearization” used by Adams [T] and Liberti m and 
also different from what is discussed in [HirnifTH] . 

In this paper, we show that the properties of QSP established in m are 
incorrect and that the algorithm they proposed need not produce an optimal 
solution. Gupta and Saxena m extended the results of m to the quadratic 
set packing and partitioning problems. These extensions also suffer from the 
same drawbacks as that of m and the algorithm in [10] could also produce 
a non-optimal solution, contrary to what is claimed. Since the algorithm of 
m is not guaranteed to produce an optimal solution, it will be interesting to 
examine its value as a heuristic. Our experimental analysis discloses that the 
algorithm of m produces good solutions for some classes of problems while 
it produces very poor solutions for other classes. 


2 The quadratic set covering problem 

Let I = {1,2,..., m} be a finite set and P = (Pi, P 2 , ■ • ■, Pn} be a family of 
subsets of I. The index set for the elements of P is denoted byJ = {1,2,...,n}. 
For each element j £ J, a cost Cj is prescribed and for each element {i,j) £ 
J X J, a cost dij is also prescribed. We refer to Cj the linear cost of the set 
Pj and c = (ci,... ,c„) the linear cost vector. Similarly dij is referred to as 
the quadratic cost corresponding to the ordered pair {Pi,Pj) and the matrix 
D = {dij)nxn is referred to as the quadratic cost matrix. 


A subset P of J is said to be a cover of /, if UjgyP^ = I. Then the 
linear set covering problem (LSP) is to find a cover L = {7r(l),..., 7r(Z)} 
such that X]i=i '^ 7 r(i) is minimized. Likewise the quadratic set covering prob¬ 
lem (QSP) is to select a cover L = {(t(1), ..., ct(Z)} such that X]i=i Ccr(i) + 
Si=l Sj=l drj(i)a(j) is minimized. 

For each i £ I, consider the vector = [an^aa ,..., Om) where 


atj 


1 i £ Pj 
0 otherwise. 
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and A = (aij)mxn be an m x n matrix. Also, consider the decision variables 
xi, X 2 , ■ ■ ■, Xn where 


1 if set Fj is selected 
0 otherwise. 


The vector of decision variables is represented by a; = {xi,... ,Xn)'^ and i 
is a column vector of size m where all entries are equal to 1. Then the LSP 
and QSP can be formulated respectively as 0-1 integer programs 


LSP: Minimize ca: 

Subject to Ax > i (1) 

X€{0,ir ( 2 ) 

and 

QSP: Minimize cx + x^Dx 

Subject to Ax > i (3) 

aie{0,ir (4) 

As indicated in m the continuous relaxation of QSP, denoted by QSP , 
is obtained by replacing the constraint x G {O,!}" by x > o, where o is the 
zero vector of size n. The family of feasible solutions of both LSP and QSP is 
denoted by S' = {a;|Aa; > i, a: G {0,1}"}. 

The following definitions are taken directly from m Any X G S is called 
a cover solution and an optimal solution to the underlying problem (LSP or 
QSP) is called an optimal cover solution. Note that each cover solution corre¬ 
sponds to a cover and vice versa. A cover V is said to be redundant if P — {j } 
for j G P is also a cover. A cover which is not redundant is called a prime 
cover. The incidence vector x that corresponds to a prime cover is called a 
prime cover solution. 

Garfinkel and Nemhauser [5] proved that if the objective function in LSP 
has a finite optimal value then there exists a prime cover solution for which 
this value is attained whenever c > o. 

Saxena and Arora claimed an extension of this result to QSP , assuming 
c > o and D is symmetric and positive semi-definite. More precisely, they 
claimed: 

Theorem 1 (Theorem 3 of m) If the ohjeetive function in QSP has fi¬ 
nite optimal value then there exists a prime cover solution where this value is 
attained. 
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This result however is not true as indicated by the following example. Let 

c=(0,0,0),A= and £>= 1^-1 1 oj 

Note that D is a symmetric and positive semi-definite matrix. For the QSP and 
QSP with A, D and c defined as above, it can be verified that x* = (1,1,1)^ 
is an optimal solution with the objective function value zero for both problems. 
The optimal cover corresponding to x* is V* = {1, 2, 3} which is a redundant 
eover since V* — {2} = {1,3} is also a cover. All other cover solutions and 
their respective objective function values are listed below: 


x^ = (1,0,1)^ redundant cover solution f(x^) = 1 

x^ = (1,1, 0)^ redundant cover solution f(x^) = 1 

x^ = (0,1,1)^ prime cover solution f(^^) = 2 

= (1,0, 0)^ prime cover solution = 2 


None of these corresponds to an optimal solution for QSP or QSP . In 
particular, no prime cover solution is optimal for the instances of QSP and 
QSP constructed above, contradicting Theorem [T] This example also shows 
that Theorem [T] cannot be corrected by replacing QSP with QSP in the the¬ 
orem. 

We now show that a variation of Theorem [T] is true, which relaxes the re¬ 
quirement of D being positive semi-definite while sign restrictions are imposed 
on its elements. This is summarized in our next theorem. 

Theorem 2 There always exists a prime cover optimal solution for QSP if c 
and D are non-negative. 

Proof Let x° G 5 be an optimal solution of QSP. Then the corresponding 
optimal objective function value is 


/(x°) = cx° + x°'"Dx° 

Let Jo be the cover corresponding to the solution x°. If Jo is a prime cover 
then statement of the theorem is correct. Otherwise we can construct a prime 
cover, let say Ji, from Jo by dropping the redundant columns. Let x^ be the 
solution of QSP with respect to the prime cover Ji and 

f(x^)=cx^-t-x^ Dx^. 

Since c and D are non-negative and Ji C Jo, 


/(*°) > /(*') 
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Since is an optimal solution to QSP, f{x^) = f{x^) and the and the 
result follows. 

The family of feasible solutions for continuous relaxations of LSP and QSP 
is represented by S' = {x\Ax > i,x > o}. The continuous relaxation of LSP 
is denoted by LSP . 

Saxena and Arora m also proposed an algorithm to solve QSP and claimed 
that it will produce an optimal solution. Their algorithm is re-stated here. 


The SaLxena-Arora algorithm for QSP 
Step 1: From the QSP, construct the corresponding QSP 
Step 2: Choose a feasible solution x'^ G S such that V/(a:°) ^ o and form the 
corresponding linear programming problem LSP as 

LSP Minimize fix^)'^x. (5) 

On solving {LSP ), let x^, be its optimal solution. Let S^ = {x^}. 

Step 3: Starting with the point x^, form the corresponding LSP , and let its 
optimal solution be x^ ^ x^. Update S^ i.e. S^ = {x^,x^}. 

Step 4: Repeat Step 3 for the point x^, and suppose at the stage S^ = 
{x^, x "^,..., a;*}. Stop, if at the {i + 1)*^ stage G S^ , then is 

the optimal solution of QSP . 

Step 5: If is an optimal solution of the form 0 or 1 then it is a solution of 
QSP otherwise, go to Step 6. 

Step 6: Apply Gomory cuts to find a solution of the 0 or 1 form and the corre¬ 
sponding prime cover. 

The algorithm discussed above suffers from various drawbacks as listed 
below. 

1. Even if c > o, and D is symmetric and positive semi-definite, the LSP 
in Step 2 could be unbounded and hence it need not have an optimal 
solution for all instances. 

2. Suppose that we apply the algorithm only for instances where LSP in 
Step 2 is bounded in all iterations. Even then, the solution produced in 
Step 4 could be non-optimal to QSP . 

3. If the algorithm terminates in Step 5 the resulting solution could be non- 
optimal to QSP. 

4. If the algorithm successfully moves to Step 6, then also the solution pro¬ 
duced could be non-optimal. 

We now illustrate each of the drawbacks discussed above using counterex¬ 
amples. 

1. Since D is a positive semi-definite matrix and c > o, the objective function 
value of QSP is bounded below by zero. However, LSP in Step 2 or in 
Step 3 need not be bounded below. Let 



6 


P. Pandey and A. P. Punnen 



/lO -3 -4 -4\ 
-3211 
-4131 


\-4 113/ 


Note that D is symmetric and positive semi-definite. Consider the instance 
of QSP with the above values for c, D, and A. Starting with the feasible 


solution x° = ( 1 , 0, 0,0)^ G S' of QSP , we get the LSP in Step 2 as 


Minimize V/(a;°)^a; = 20a;i — 6 x 2 — 8 x 3 — 8 x 4 
Subject to: xi + X 2 > 1 


Xi + X 3 > 1 
Xi + X 4 > 1 
Xj > 0 for j = 1 , 2 ,3,4. 


This problem is unbounded. Thus the algorithm can not be applied in this 
case. The immediate conclusion is that the Saxena-Arora m algorithm 
is potentially applicable only to those QSP instances where the resulting 
LSP is bounded in every step. 

2. The algorithm can fail in Step 4. The Saxena-Arora algorithm claims to 
produce an optimal solution of QSP in Step 4 but this may not be true 
always. Consider the data 




c = (0,0,0, 0), A = 10 10 and D = 


Note that D is symmetric and positive semi-definite. Consider the instance 
of QSP with the above values for c, D, and A, we get the QSP as 

Min f{x) = 10xl + 8 x 2 + 8 x\ + 'ix\ + ‘ixiX 2 +^xiX 3 + ‘ixiX 4 -I- 2 X 2 X 3 + 2 X 2 X 4 + 2 X 3 X 4 
st: xi + X 2 > 1 
Xi + X 3 > 1 
Xi + X 4 > 1 


Xj G {0,1} for j = 1, 2, 3,4. 


and 


Vf{x) =( 20 x 1 + 4 x 2 + 4 x 3 + 4 x 4 , 6 x 2 -I- 4xi -I- 2 x 3 -I- 2 x 4 , 
6 x 3 -I- 4xi -I- 2 x 2 + 2 x 4 , 8 x 4 -b 4xi -b 2 x 2 + 2 x 3 )^ 
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Select the feasible solution = (1,0,0,0)^ G S' of QSP . Construct the 
LSP with respect to a;°, the objective function of LSP is 


Vf{x'^)^x =20xi + 4x2 + 4x3 + 4x4 

Note that x^ = (0,1,1,!)^ is an optimal solution to this LSP . Thus, 
we set S^ = {x^}. Now, using x^ construct the new LSP , and the opti¬ 
mal solution to this LSP is x^ = (1,0,0, 0)^. Since x^ ^ S^, construct 
the new LSP , the optimal solution to this LSP is x^ = (0,1,1,!)^. 
Since x^ G S^, in Step 4 the algorithm concludes that x^ is an opti¬ 
mal solution of QSP with objective function value 16. However, x* = 
(0.714286, 0.285714,0.285714,0.285714)'^ which is a better solution for QSP ', 
contradicting the optimality of x^. Thus the algorithm could fail in Step 4. 

3. As per the Saxena-Arora algorithm. Step 5 produces an optimal solution 
to QSP and if this optimal solution is binary, they claim this solution 
to be an optimal solution of QSP. We now show that a binary solution 
produced in Step 5 need not be an optimal solution to QSP. For example. 

Consider the data 


c= (0,0,0,0), A 


/I 1 0 0\ 
10 10 
\1 0 0 1 / 


and 


D = 


/4 1 1 1\ 
12 0 0 
10 2 0 
\1 0 0 2 / 


Note that ZD is a symmetric, positive semi-definite and non-negative. As 
noted in Theorem [51 a prime cover optimal solution exists for this QSP. 
But still the Saxena-Arora algorithm fails to produce an optimal solution 
for QSP. Consider the instance of QSP with the above values for c, ZD, and 

A. 


Select the feasible solution a;° = (|, 4, 4, 4)'^ G S which is also an optimal 
solution of QSP . Construct the LSP with respect to x^ and the objective 
function is Vf{x^)'^x = ^Xi -|- |x2 -|- |x3 -I- |x4. 

x^ = (0,1,1,1)^ is an optimal solution to this LSP . Thus, we set S^ = 
{x^}. Now, using x^, construct the new LSP with the objective func¬ 
tion as V/(x^)^x = 6xi -I- 4x2 + 4x3 + 4x4. An optimal solution to this 
LSP is x^ = (1,0,0, 0)^. Since x^ ^ S^, we update S^ = {x^,x^}. Start¬ 
ing with x^, construct the LSP with the objective function V/(x^)^x = 
8 xi- 1 - 2 x 2 - 1 - 2 x 3 - 1 - 2 x 4 . An optimal solution to this LSP is x^ = (0,1,1, l)"^. 
Since x^ G 5”^, the algorithm concludes that x^ is an optimal solution of 
QSP . Since x^ contains 0 and 1 entries only, as per the algorithm, it is an 
optimal solution to QSP and the corresponding objective function value is 
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However x* = (1,0, 0,0)^ is a better solution to the QSP with objective 
function value f{x*) = 4. Thus, the solution produced by the the Saxena- 
Arora algorithm for the above instance of QSP is not optimal. 

In the previous example if = (0,1,1,1)^ is selected instead of 
the algorithm produces xi = (1, 0, 0,0), X 2 = (0,1,1,1), and xs = (1,0, 0,0), 
leading to an accurate optimal solution xs = (1,0,0, 0) to QSP. Note that 
a;o = (0,1,1,!)^ and (1,0, 0,0)^ are alternate optimal solutions of LSP 
with the objective function Vf{x^)'^x = ^xi + ^X 2 + ^xs + 5 X 4 . It is 
easy to show that trouble of the Saxena-Arora algorithm is not because of 
the presence of alternate optimal solutions, leading to a choice in selection. 
This can be demonstrated with the same example but by selecting a dif¬ 
ferent starting point as given below. 

Select the feasible solution = (l,i,0,0)^ € S' of QSP . Construct 
the LSP with respect to x^ and the objective function is Vf(x^)"^x = 
9xi -I- 4 x 2 + 2 x 3 + 2 x 4 . x^ = (0,1,1,1)^ is the unique optimal solution to 
this LSP (easily verifiable by enumerating the basic feasible solutions). 
Thus, we set S^ = {x^}. Now, using x^, construct the new LSP with the 
objective function as Vf(x^)^x = 6 x 1 -I- 4 x 2 + 4 x 3 + 4 x 4 . The unique op¬ 
timal solution to this LSP is x^ = (1,0,0,0)^. Since x^ ^ S^, we update 
S^ = {x^,x'^}. Starting with construct the LSP with the objective 
function Vf{x‘^)'^x = 8 x 1 -I- 2 x 2 + 2 x 3 -I- 2 x 4 . The unique optimal solution 
to this LSP is x^ = (0,1,1,1)^. Since x^ G S^, the algorithm concludes 
that x^ is an optimal solution of QSP . Since x^ contains 0 and 1 entries 
only, as per the algorithm, it is an optimal solution to QSP and the corre¬ 
sponding objective function value is 6 . The solution produced by the the 
Saxena-Arora algorithm for the above instance of QSP is not optimal. 

4. As per the Saxena-Arora algorithm. Step 5 produces an optimal solution 
to QSP and if this optimal solution is not binary, the algorithm proceeds 
to Step 6 where Gomory cuts are applied to find a solution which they 
claim to be an optimal solution to QSP. We now show that Step 6 need 
not produce an optimal solution to QSP even if the solution produced in 
Step 4 is optimal for QSP . 

In point 3 we gave a counterexample where the solution is a basic feasible 
solution (BPS) to LSP which is binary but not optimal to QSP . Note 
that QSP is a continuous quadratic problem and an optimal solution need 
not correspond to an extreme point. We now illustrate that if the LSP 
solver works with any solution (such as interior point methods) and not 
necessarily with BPS (as in simplex method) it may be possible to get an 
optimal solution to QSP in Step 4. Por example 
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Consider the instance of QSP from the previous case. = (f, I, S 

S is an optimal solution of QSP . Construct the LSP with respect to 
and the resulting objective function is Vf{x°)'^x = ^xi + |x 2 + fxa + |x 4 . 
The algorithm produces xi = , X 2 = leading to 

an accurate optimal solution X 2 = (|, \) to QSP and the algorithm 

successfully completes Step 4. 

To apply Gomory cut,first reduce the non-basic feasible solution (non-BFS) 
to a basic feasible solution (BPS). From the previous exampl, the optimal 
non-BFS (|, \) of LSP to an optimal BFS x^ = (0,1,1,1)^ of LSP . 

Since this is binary , no cutting plane will be added and a Gomory cut phase 
terminates with the non-optimal solution x^ = (0,1,1,1)^ of QSP. 

Alternatively if we do not reduce the non-BFS to a BFS to apply Gomory 
cuts, but use any Integer programming (IP) solver to compute an optimal 
integer solution to LSP we could still get non-optimal solution. For ex¬ 
ample: Solving the LSP at (|, j, j) for 0-1 optimal solution we could 
get x^ = (0,1,1,1)^ as the optimal 0-1 solution of LSP . This is not an 
optimal solution to QSP. (We note that the paper [T7] does not say any¬ 
thing about the use of general IP solver; but we mentioned it here for the 
clarity and completeness). 


3 The quadratic set packing and partitioning problems 

A subset FI of J is said to be a pack of I if Ujg// Pj = for j, k € H, j ^ 

k, implies Pj Q\Pk = ^- Then the linear set packing problem (LSPP) is to select 
a pack V = {7r(l),..., 7r(u)} such that maximized. Likewise, the 

quadratic set packing problem (QSPP) is to select a pack L = {cr(l),..., cr{l)} 
such that c^(d + ELi T,j=i i)a(j) is maximized. 

Let A = {aij)mxn be as defined in Section [2] Also, consider the decision 
variables xi,X 2 , ■ ■ ■ ,Xn where 


1 if j is in the pack 
0 otherwise. 


The vector of decision variables is represented as a; = (xi,..., x„)^ . Then 
the LSPP and QSPP can be formulated respectively as 0-1 integer programs 


LSPP: Maximize ca: 


Subject to Ax < i 

a; G {0,ir 


( 6 ) 

(7) 
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and 


QSPP: Maximize cx + x^Dx 

Subject to Ax < i 

a: e {0,ir 


( 8 ) 

(9) 


The continuous relaxations of LSPP and QSPP, denoted respectively by 
LSPP(C) and QSPP(C), are obtained by replacing the constraint x S {0,1}" 
by a; > o, respectively in LSPP and QSPP. 

The family of feasible solutions of both LSPP and QSPP is denoted by 
S = {x\Ax < i,a; S {0,1}"} and the family of feasible solutions for their 
continuous relaxations is denoted hy S = {x\Ax < i, a; > o}. 

Following are some definitions given in m- A solution X € S which satis¬ 
fies (HJ and (HI) is said to be a pack solution. For any pack V, a column of A 
corresponding to j* G F is said to be redundant if F — {}*} is also a pack. If 
a pack corresponds to one or more redundant columns, it is called a redundant 
pack. A pack V* is said to be a prime pack, if none of the columns correspond¬ 
ing to j* G F* is redundant. A solution corresponding to the prime pack is 
called a prime packing solution. 

From the definition of a redundant column given above (as in nni) , zero 
vector is the only prime packing solution for the set packing problem. Thus 
the results of m are incorrect with respect to their definitions. We believe 
the sign in the above definition of redundant column discussed in |10j is 
a typo and it is probably supposed to be “ U ” which is consistent with the 
definitions given in m by the same authors. Hereafter, we use this modihed 
definition. 

Thus, for any pack V, a column of A corresponding to j G J is said to be 
redundant if F U {}} is also a pack. If a pack contains one or more redundant 
columns, it is called a redundant pack. A pack V* is said to be a prime pack, 
if none of the columns corresponding to j G J is redundant. A solution corre¬ 
sponding to the prime pack is called a prime packing solution. 

Gupta and Saxena [10] assumed Z) to be a negative semi-definite matrix 
and extended most of the results for QSP in m to QSPP. In particular, they 
claimed that: 

Theorem 3 (Theorem 2 of [10]) If the objective function in QSPP has 
finite value then there exists a prime packing solution where this value is at¬ 
tained. 

Because of the definition of the prime pack solution given by Gupta and 
Saxena m, a prime pack is always a zero vector hence the theorem is given 
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incorrect. The theorem is still incorrect even if we use the modified definition 
m which is indicated above. 

For example, consider an instance of QSPP with 


c= (0,0,0), A = 

Note that D is symmetric and negative semi-definite. 




■-2 1 O' 

'1 1 o' 

,D = 

1 -2 1 

10 1 



0 1-2 


X* = (0, 0, 0)^ is an optimal solution for the QSPP with the objective 
function value zero. We list below all prime pack solutions with the objective 
function values. 

— {1,0, 0} prime pack solution = —2 

x'^ = {0,1,1} prime pack solution f{x^) = —2 
Note that none of these solutions are optimal. 

We now show that a variation of Theorem [3] is true and this is summarized 
in our next theorem. 


Theorem 4 There always exists a prime pack optimal solution for QSPP if 
all elements of c and D are non-negative. 

Proof Let a;° G be an optimal solution of QSPP. Then the corresponding 
optimal objective function value is 


f{x°) =cx° + x°'"Dx° 

Let Jo be the pack corresponding to the solution x°. If Jo is a prime pack 
then we are done. Otherwise we can construct a prime pack, let say Ji, from Jo 
by adding the redundant columns. Let x^ be the solution of QSP with respect 
to the prime pack Ji and 


f {x^) = cx^-\-x^ Dx^. 

Since Ji obtained by adding redundant columns to Jq, therefore, Jq C Ji. 
When all elements of c and D are non-negative, or 


f{x°) < f{x^) 


Since x^ is an optimal solution to QSPP, f{x^) = f{x^) and the proof 
follows. 
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Along the same lines as in m, the authors of m, provide a solution al¬ 
gorithm for QSPP. Following the insight generated in our counter examples in 
Section [21 and by the above observation, it is not difficult to construct counter 
examples to show that the algorithm of nni need not provide an optimal so¬ 
lution for QSPP. 

If in equation ([5]) we replace constraints Ax < i with Ax = i, then 
QSPP changes into quadratic set partitioning problem. Gupta and Saxena 
[10) proposed a similar algorithm for the quadratic set partitioning problem, 
which has similar issues as in the quadratic set packing problem. We omit the 
discussion about the quadratic set partitioning problem. 


4 Computational results 

Since the algorithm of m is not guaranteed to be optimal, it would be in¬ 
teresting to examine its value as a heuristic to solve QSP. We have conducted 
some preliminary experimental analysis to assess the value of the Saxena-Arora 
algorithm as a heuristic using different classes of test problems. 

The test data was taken from standard benchmark problems for the set 
covering problem [211411^ . and the vertex covering problem |18j . with appro¬ 
priate amendments to incorporate quadratic objective. In this class, we took 
only small size instances since the quadratic problem is much more difficult 
and time consuming to solve compared to their linear counterparts. We have 
also generated some quadratic vertex cover instances on random graphs taken 
from m- We divided computational experiments into two different categories, 
with each c > 0, while in category 1: £) is a positive semi-definite matrix and 
in category 2: D is non-negative and positive semi-definite. 

Each element of the linear cost vector c is a random integer from the in¬ 
terval [3,5]. Since the quadratic cost matrix D is positive semi-definite, there 
exists a square matrix B such that D ~ . This D is generated by a ran¬ 

dom square matrix B where each element of B is a random integer between 
-10 and 10. When D is non-negative and positive semi-definite, each element 
of B is selected as a random integer between 0 and 20. 

The Saxena-Arora algorithm was coded in C-I--I- and tested on a PC with 
windows 7 operating system, Intel 3770 i7 3.40 GHz processor and with 16 
GB of RAM. We also used CPLEX 0-1 integer quadratic solver (version 12.5) 
to compute exact (heuristic) solutions. For each instance that we tested, we 
set CPLEX time limit to be the same as the time taken by Saxena-Arora 
algorithm and also run CPLEX by doubling this running time. These two im¬ 
plementations provide heuristic solutions and were compared with the solution 
produced by the Saxena-Arora algorithm. 
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In the tables, ti is the cpu time taken by Saxena-Arora algorithm. The col¬ 
umn “CPLEX Sol (ti)” represents the heuristic solution obtained by CPLEX 
by fixing its running time to ti and the column “CPLEX Sol (2ti)” represents 
CPLEX run with 2ti upper bound on the execution time. The column “neg¬ 
ative entries in Q” provides percentage of negative entries in the matrix D. 
The column “Sol” refers the objective function values. CPLEX quadratic solver 
takes more time to solve QSP to optimality when D is positive semi-definite 
compare to the instances when D is non-negative and positive semi-definite. 
Therefore, Table [1] reports lower bound value and Table [5] reports optimal 
solution value. 


Table 1 Benchmark instances, D is positive semi-definite 


problem 

size 

m 

n 

Lower 

bound on 
opt 

Saxena Algo. 

cpu time ti g ^ 

(sec) 

CPLEX Sol 
(ii) 

CPLEX Sol 

(2ti) 

Negative 
entries in 
D(%) 

QscpcvcOG 

192 

240 

48 

156 

241189 

71 

70 

43.59 

Qscpcyc07 

448 

672 

112 

125 

1486218 

192 

187 

46.00 

Qscp41 

1000 

200 

432.27 

78 

7743176 

455 

455 

44.48 

QscpeS 

500 

50 

3.307 

63 

1777128 

10 

10 

46.46 

Qscpe4 

500 

50 

3.455 

62 

1638640 

10 

10 

45.69 

Qscpe5 

500 

50 

3.393 

47 

1918480 

28 

15 

45.69 

O.Kraph50-01 

612 

50 

111.75 

156 

144 

122 

122 

41.04 

O.Kraph50-02 

490 

50 

93.6684 

172 

641 

124 

124 

42 

O.Kraph50-03 

735 

50 

101.25 

156 

153 

152 

152 

40.23 

O.eraph50-04 

612 

50 

87.9275 

156 

1057 

149 

149 

45.6 

O.Kraph50-05 

490 

50 

108.25 

156 

127 

124 

124 

45.6 

O.eraph50-06 

857 

50 

89.71486 

187 

261 

154 

154 

41.6 

O.eraph50-07 

735 

50 

85.1275 

156 

1738 

119 

119 

48 

O.eraph50-08 

612 

50 

89.1389 

156 

1030 

100 

100 

44 

O.eraph50-09 

980 

50 

89.6389 

156 

600 

146 

146 

47.52 

O.eraph50-10 

612 

50 

92.5273 

156 

2363 

142 

142 

48.72 

Qfrb30-15-1 

17827 

450 

810.306 

11185 

5106 

1494 

1494 

42.56 

Qfrb30-15-2 

17874 

450 

799.139 

7909 

18906 

1479 

1479 

44.59 

Qfrb30-15-3 

17809 

450 

799.139 

6100 

11497 

1478 

1477 

44.59 

Qfrb30-15-4 

17831 

450 

787.806 

5274 

21930 

1475 

1475 

44.46 

Qfrb30-15-5 

17794 

450 

799.139 

5616 

14032 

1478 

1475 

44.59 


When £> is a random positive semi-definite matrix. Table [T] shows that 
the Saxena-Arora algorithm does not return a good quality solution for QSP. 
Note that a general purpose solver like CPLEX obtained much better solu¬ 
tions within the same time limit for the test problems used. But when D is 
non-negative and positive semi-definite. Table [5] shows that the Saxena-Arora 
algorithm produced solutions as good as those produced by CPLEX for many 
instances. For vertex cover instances the Saxena-Arora algorithm produced an 
optimal solution. For the set cover instances CPLEX produces better solutions 
than the Saxena-Arora algorithm. Thus, for D is non-negative and positive 
semi-definite, the Saxena-Arora algorithm could be used as a heuristic to solve 
QSP. As our counter example indicates, even for this class the Saxena-Arora 
algorithm need not produce an optimal solution. 
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Table 2 Benchmark instances, D is non 

-negative and positive semi-definite 


problem 

size 

Optimal sol 

Saxena Algo. 

CPLEX Sol 
(ti) 

CPLEX Sol 

(2ti) 

m 

n 

cpu time ti 
sec 

Sol 

Qscpcyc06 

192 

240 

147523 

483 

147523 

147523 

147523 

Qscpcyc07 

448 

672 

726070 

780 

726070 

726070 

726070 

Qscp41 

1000 

200 

7271 

343 

7271 

7271 

7271 

QscpeS 

500 

50 

7 

390 

3369 

7 

7 

Qscpe4 

500 

50 

8 

515 

3726 

8 

8 

Qscpe5 

500 

50 

7 

359 

1854 

7 

7 

Qgraph50-01 

612 

50 

5149 

187 

5149 

5149 

5149 

Qgraph50-02 

490 

50 

5108 

156 

5108 

5108 

5108 

Qgraph50-03 

735 

50 

5163 

188 

5163 

5163 

5163 

Qgraph50-04 

612 

50 

9272 

156 

9272 

9272 

9272 

Qgraph50-05 

490 

50 

4168 

125 

4168 

4168 

4168 

Qgraph50-06 

857 

50 

8042 

124 

8042 

8042 

8042 

Qgraph50-07 

735 

50 

7088 

124 

7088 

7088 

7088 

Qgraph50-08 

612 

50 

4670 

109 

4670 

4670 

4670 

Qgraph50-09 

980 

50 

8872 

141 

8872 

8872 

8872 

Qgraph50-10 

612 

50 

6614 

125 

6614 

6614 

6614 
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